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We investigate the evolutionary aspects of some integrable soliton models whose La- 
grangians are derived from the pullback of a volume-form to a two-dimensional target space. 
These models are known to have infinitely many conserved quantities and support various 
types of exact analytic solutions with nontrivial topology. In particular, we show that, in spite 
of the fact that they admit nice smooth solutions, wave propagation about these solutions 
will always be ill-posed. This is related to the fact that the corresponding Euler-Lagrange 
equations are not of hyperbolic type. 


I. INTRODUCTION 


Nonlinear theories supporting topologically 
protected solitons appear in a variety of phys¬ 
ical situations, from nuclear interactions to con- 

dense. 

with a source space-time (M, g) and consider 
fields if taking values not in a vector space, but 
rather in a smooth target Riemanniann manifold 
(N,/i). The topological character of solutions 
appears explicitly when we look for static con¬ 
figurations of finite self-energy. Then, hypothe¬ 
sis concerning the behaviour of fields at spatial 
infinity effectively imply that they can be classi¬ 
fied according to some homotopy class. Roughly, 
this means that solutions are disjoint and can¬ 
not be smoothly deformed into each other due to 


topological reasons. Relevant examples of soli¬ 


ton t 


els [ 21 , the SU{2) Skyrme model of pions |J], and 
the Faddeev-Niemi model of knotted hopfions 


leories of this type include the 0{k) a mod- 


s, 


B 


The point of this letter is to shed light into 
the evolutionary aspects of some integrable soli¬ 
ton models whose Lagrangians are derived from 
the pullback ip*e of a volume-form e to a tar¬ 
get space (see, for instance, B)- These mod¬ 
els are important for various reasons: i) They 
generalize to higher space-time dimensions some 
concepts which proved useful in the context of 
(l-l-l)-dimensional field theories such as the zero 

n 

curvature of Lax-Zakharov-Shabat [7[. ii) The 
corresponding system of quasilinear second order 
partial differential equations (PDE’s) B is some- 
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what simpler and often enables analytical treat¬ 
ment of soliton dynamics, hi) They can play a 
role in the phenomenological description of phys¬ 
ically relevant systems such as non-perturbative 
regimes of QCD and BPS models of baryons at 
large Nc js'lll- Besides - and perhaps even more 
persuasively - the models are known to have in¬ 
finitely many conserved quantities and support 
various types of exact solutions with nontrivial 
topology. In particular, they predict the exis¬ 
tence of static and time-dependent hopfions in 

(1-1-3) space-time dimensions and exotic tex- 

□ 

tures for higher dimensions |12l |. 

The main question we address here is whether 
the notion of integrability (in the generalised 


sense of 


Q] and Q]) i 


is compatible with the hyper¬ 


bolic nature one usually expects in a relativistic 
field theory [l^ . For the sake of simplicity we re¬ 
strict ourselves to the case of a two-dimensional 
target space. In particular, we show that, de¬ 
spite the equations of motion admitting infinite 
exact analytical solutions, their linearised ver¬ 
sions are always ill-behaved. In other words, 
small disturbances around a given solution will, 
in general, propagate in an undesirable way, sig¬ 
nalling to instabilities. We shall see that this 
result does not depend on the topology of the 
target manifold but, rather, is due to the spe- 
cihc type of nonlinearities present in the Euler- 
Lagrange equations. In particular, a detailed 
analysis reveals that the associated PDF’s are 
not of hyperbolic type. 

The organization of this paper is as follows. 


In section II we review some geometrical aspects 
of the Lagrangian theory of maps and establish 
the basic mathematical machinery. Section III 
deals with maps whose Lagrangians are derived 
from the pullback of a volume-form to a target 
space. Here we start with a general framework 
and then restrict ourselves to the simpler case 
of a two-dimensional target space. We continue 
in section IV by investigating the characteris¬ 
tic surfaces of the model and discussing the evo¬ 
lutionary aspects of the linearized equations of 
motion. We conclude in section V by mention¬ 
ing a number of possible extensions that could 
lead to a well behaved theory. 


II. GEOMETRICAL FRAMEWORK 


We start, briefly, placing the theory in the 
general geometrical setting of harmonic/wave 
mappings. Although this construction is hardly 
employed in the context of solitons, we believe 
that it is particularly adequate for our purposes 
and sets the background for future related prob¬ 
lems. Our steps here proceed very much in 
the same way as Fells and Sampson [l6| and 

n n 

Choquet-Bruhat [171] (see also [18| for some phys¬ 
ical motivations). 


A. Kinematics 

Write g) for a (1 -|- m)-dimensional 

space-time with metric signature (-I-, —, —,...), 
(N”, h) for a n-dimensional riemannian mani- 





B. Dynamics 
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fold with metric h and consider smooth maps of 
the form 

(f : ^ N". (1) 

If (o = 0, and {A = de¬ 

note local coordinates in the base and target, 
respectively, the map reads The 

differential of (/9 at x G is, in a sense, the 

best linear approximation of (p near x, 


When m > n, as is always the case for the 
integrable theories we shall discuss in the next 
sections, it follows 0 < rank(L) < n and cjfc = 0 
for all k > n. In this situation the most general 
first order action is provided by 

5[v?] = / an,(p) dvg, (6) 

J^l + m 


dp,, : ^ T^(x)N", (2) 

and can be used to pullback tensors living in N” 
to In particular, using the convention 

daP>^ = the pulled back metric pi*h reads 

Lab{x) = hAB{p>{x))p>\p>\ (3) 


with dvg the element of volume in The 

Euler-Lagrange equations can be written in a 
compact form if we introduce a linear connection 
in the associated vector bundle E = (g) 

If u\{x) is a smooth cross-section of 
E, its covariant derivative is written as |18l |. 


The eigenvalues of (which is often called the 


strain for the map 


201) at X G solve the 


characteristic polynomial 


l+m 

det(L“, - X5\) = Y, = 0, (4) 

k=0 

where is the fe-th elementary symmetric func¬ 
tion of L“. 


n ^ (5) 

with the generalized Kronecker tensor 

and (To = 0. The basic idea is to use the in¬ 
variants (Tfe as ingredients for a general Lorentz 
invariant Lagrangian prescription. In particular, 
defining L := 


= dbu\ - 

(7) 

where and are the Christof- 

fel symbols corresponding to the the metrics g 
and h, respectively. Generalization to mixed ob¬ 
jects with more than one internal index (A) fol¬ 
lows the same route: simply add one connection 
{N)pA^^ term to ‘covariantise’ each of the inter¬ 


nal indices 


Q- 


With this notation the critical points of 
satisfy the nonlinear equations of motion 

dC 




2 d(p 


B ’ 


( 8 ) 


‘Xi — [L], (72 — ([L]^ — [L^])/2, where 

a3 = ([Lf-3[L][L2] + 2[L3])/6, 

with the brackets [ ] denoting trace operation, 
for conciseness. 





.b bi...bk-i 
Oa ai...afc_i 


Ij^i 

bl'" 6fe-i 




(9) 
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and Ck = dCjdak for conciseness. The pattern 
here is simple: the k-th term inside brackets in¬ 
volves (fc — 1) powers of the strain. Written more 
explicitly, Eq. (I8|) becomes 

uAB ptn 


i.e. C C{H) with H = In¬ 

terestingly, the associated action is a particular 
instance of ([6]) as it directly follows from defini¬ 
tions dH) and (fTOl) that H = nlu^. 

Thus, for theories living in a space-time 
with m > n, one considers the family 
of Lagrangians 


revealing that various types of nonlinearities can 
be present. Generically, the first term leads to 
quasi-linear contributions while the others give 
rise to semi-linear terms. In particular, the r.h.s. 
plays the role of a potential and can be cho¬ 
sen according to different motivations. Note also 
that the PDE’s are covariant with respect to co¬ 
ordinate re-parametrizations both in and 

N". 

III. PULLBACK OF THE VOLUME 

FORM AND INTEGRABLE MODELS 

A. n-dimensional targets 

A class of nonlinear theories which well fits 
in the geometrical formalism described so far is 
one where the Lagrangian consists of a function 
of the pullback ip*e of the pertinent volume form 
on a given target N”, 

( 10 ) 

with eAi...An representing the totally antisym¬ 
metric Levi-Civita tensor of dimension n. Specif¬ 
ically, one considers a Lagrangian that is propor¬ 
tional to a smooth function of the square of ip*e 


S[.p\int = / ^(cTn) dVg. (11) 


The restriction m > n is important since in 
the case where the number of spatial dimensions 
m coincides with the dimension of the target 
n the corresponding static equations are triv¬ 
ial ^ . When m is further restricted to the 
form C{an) oc an, with q a positive real parame¬ 
ter, it is possible to choose q so as to avoid Der¬ 
rick’s scaling arguments and various types of ex¬ 
act solitonic solutions were obtained in the case 
of spherical targets = S"' (see, jb] and refer¬ 
ences therein). These solutions are topological 
in the sense that they’re characterised by the 
corresponding homotopy group 7rm(§"') G Z. 

Models of this form are integrable in the sense 
that they posses a generalised zero-curvature 
representation and an infinite number of local 
conservation laws discussed in a series of papers 
0i) 0j) Qj Ql- Roughly, the Noether currents 
are associated to the huge amount of symmetry 
implied by the invariance of the theory under 
volume preserving diffeomorphisms on the target 
space. Using ([9]) we obtain, after some algebra 


iLl oc Cne^ 


TT ai...an,.Ai An-\ 

nt Ai...An-i^a ^ a„_i • 
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Note that, in this form, both the volume-form 
e and its pullback ip*e appear explicitly. Now, 
from 0 it directly follows 

= 0 D[,56](^^ = 0, (12) 

which imply that the equations of motion for the 
integrable models read 

(13) 

where, now, V“ denotes ordinary covariant dif¬ 
ferentiation in M. 


B. 2-dimensional targets 


From now on we shall focus on integrable 
models living in a (3-|-l)-dimensional space-time 
and taking values in a 2-dimensional target 
space. Apart from these restrictions, we make 
no further hypothesis on the geometry/topology 
of and N^. One first construction of this 

type was considered by Aratyn, Ferreira and 
Zimerman (AFZ) with base space and tar¬ 
get the symmetric space SU{2)/U{1) = S^. As¬ 
suming the fractional Lagrangian 




a. 


3/4 


dv„ 


/R1+3 


(14) 


so as to bypass the usual obstacle due to Der¬ 
rick’s scalings arguments, they show that Eq. 
m is solvable and infinitely many static analyt¬ 
ical solutions of toroidal symmetry, classifie d by 


a Hopf index Qh = 7r3(S^) G Z, exist 2]|, 2^. 


A second model, which explores the invariance 


under the conformal group 50(4,2), is provided 
by the strongly-coupled model 


S[(p]int = -7: 0-2 dVg, (15) 

4 jRl+3 


introduced by Ferreira in 


23] (see also 24l]). 


In this case, although the solitons cannot be 
brought to rest there exist exact nontrivial time- 
dependent hopfions carrying angular momen¬ 
tum. As stressed in 2^, this is a rare example of 


an integrable theory in four dimensions and its 
solitons may have some role in the low energy 
limit of Yang-Mills theory. 

In principle, the assumption of a flat space- 
time is not mandatory and more ambitious mod¬ 
els using curved source spaces have been consid¬ 
ered. In 25(, for instance, the authors consid¬ 
ered a space-time of the form x M and the 
Poincare disk as a target. The important point 
here is that, for all of these models, Eq. m 
reduces to 


da (^/^ = 0, (16) 

with Hab = eAB^\‘P^b the area two- 

form on N^. The following three facts hold for 
Eq. ()16l) : i) it consists a two-dimensional system 
of second-order, quasi-linear PDE’s, for the map; 
ii) the connection of the target space 
does not appear explicitly in the equations, im¬ 
plying that various types of self-interactions do 
not contribute, iii) if (p^{x) is a solution of (I13p 
for a given Hab it is also a solution for Hab if 
det(/iAB) = det(hAs). 
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IV. CHARACTERISTIC ANALYSIS 


The natural question to ask for equations (fTUl) 
is whether they pose an initial value formula¬ 
tion: in other words, we would like to deter¬ 
mine the complete map ip : (N^, h) 

from its values p\y. and first derivatives dp\Y, re¬ 
stricted to a space-like sub-manifold S C M. 
For well-posedness to hold, it is crucial that 
small disturbances (or, equivalently, discontinu¬ 
ities) in the initial data propagate in a con¬ 
trolled/predictable way with a definite finite ve¬ 
locity of propagation. Roughly speaking, this 
means that the associated characteristic surfaces 
about a background solution have to be de¬ 
scribed by algebraic varieties with the topology 


of convex cones for each point of (see 261 


and 


27l | for a similar analysis in the contexts of 


the Baby-Skyrme and Faddeev models ). 

In order to clarify this point, it is convenient 
to rewrite Eq. m as 

d^) dadbP^ + JaIp, dp) = 0, (17) 


where JAip,dp) stands for semilinear terms in 
p (whose explicit form is unnecessary for our 
discussion) and so-called principal 

part of the system. Here, = 0, 

with brackets denoting antisymmetrization. A 
closer inspection of Eq. (USD gives 

^'^°‘^AB = i^AP^BQ — apNbq) d'^p^^p^^ — 

-Q^^Nab. (18) 

with Nab = eAP^BqdcP^d’^p^ and ^ = 
2C22IC.2- As it is well known, the principal 


part almost completely controls the qualitative 
behaviour of solutions of a PDE 2^. In par¬ 


ticular, it determines the evolution of linearized 
waves about some smooth background solution 
Poix) and gives rise to the causal structure of 
the theory on top of this solution. Although 
some progress has been made in this direction 
for the Skyrme and related models (see, for ex¬ 


ample, Wong’s paper 


291]), to the best of our 


knowledge no other authors have considered the 
detailed structure of (fT71) in the case of integrable 
soliton models. 

Physically, the characteristics can be iden¬ 
tified with the infinite-momentum limit of 
the eikonal approximation (or, equivalently, 
with the surfaces of discontinuity obtained via 
Hadamard’s method 3^). Briefly, they can be 


obtained as follows: recall that, for a covector 
ka G TJM, the principal symbol is given by the 
contraction (see, for instance, [31[) 


MAB{k) = M^^ABkakb. 


(19) 


In our case, we have 

MAB{k) = \^\\hAB - \k\^gNAB - ^A^B “ 

-^NApNBQe^i^ ( 20 ) 

with = d°'4>^ka and \l\\, j/cj^ the usual norms 
with respect to the corresponding metrics. Now, 
If p{){x) is a smooth solution of (fTTp . the deter¬ 
minant \MAB{k)\ is a real valued function of x 
and k and is called the characteristic polynomial 
V{x,k) around the solution po{x). The set 

N^ = {keT*M\V{x,k) =0,k^0} (21) 
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is called the characteristic set and consists of the 
locus of normals to the characteristic surfaces at 
X. As Eq. ([201) is given by a 2 x 2 matrix, the 
characteristic polynomial reduces to the form 

V{X, k) = (22) 

with a completely symmetric tensor den¬ 

sity depending on ifo{x). Thus, the wave nor¬ 
mals are determined by the vanishing sets of 
a multivariate polynomial of fourth order in ka 
in the cotangent space. The resulting algebraic 
variety changes from point to point in a way 
completely prescribed by the background solu¬ 
tion and the nonlinearities of the given integrable 
model. For such theories (as for quasi-linear the¬ 
ories in general) wave velocities are not given a 
priori, but change as functions of initial data, 
directions of propagation and wave polarization 
38|. 

It is quite common in nonlinear field theories 
described by second-order PDE’s that the char¬ 
acteristic polynomial factorizes into simpler irre¬ 
ducible polynomials. This phenomena happens, 
for instance, in the context of nonlinear electro¬ 
dynamics 32 


inside media 
ity 


Born-Infeld theories 


33l |. optics 


and Lovelock theories of grav- 


35l |. The integrable theory of maps into a 


with h = det(/iyiB), 


Viix, k) = G^f^^kah, V 2 (x, k) = G'{ 2 )kckd, 


and 


/^ab — ab , Tja zrcb 
Lr(i) = 0 - 25 ' + H , 

G ab ab /■ rja rrcb 

(2) = 5 “ S-H c-H • 


( 2 )' 


(24) 

(25) 


In order to derive these relations we used the 
fact that, for a 2 -dimensional target and a four 
dimensional base, era = CT 4 = 0. In this case 
the Cayley-Hamilton theorem then implies the 
relation = aiL^ — (T 2 L^ which is important 


for the factorization 


36|. 


At first sight it appears that the variety of 
wave normals is given by a product of quadratic 
Lorentzian cones in T*M. What about the char¬ 
acteristic surfaces themselves? The theory of 
PDE’s proceeds by showing that they are gener¬ 
ated by bi-characteristic rays a:“(A). In the case 
where V factorizes, then we must use each Vi 
(i=l, 2 ) instead of V in defining the character¬ 
istics. Consequently, the bi-characteristic rays 
{ 3 :“(A)} are solutions of the canonical equations 


m 

dka 


ka = - 


m 

dx°- 


(^ = 1 , 2 ), 


two-dimensional surface is not an exception. In¬ 
deed, for the matrix (( 201 ) a tedious, but straight¬ 
forward, calculation reveals that the character¬ 
istic equation always factorizes into a product of 
quadratic terms, i.e. 


V{x^ k) = hVi{x, k)V 2 {x, k) = 0, 


(23) 


where the dot means derivative with respect to 
the parameter A. The idea here is to use the 
first equation to solve for ka in terms of and 
use the second to obtain a second order equa¬ 
tion for x“(A), as is common in Hamiltonian me¬ 
chanics/optics. It is clear that this can only be 
achieved if both the quadratic forms G'^-^ (which 
depend on the background) are nondegenerate. 
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Unfortunately, this is not true for the ex¬ 
pressions ([2^ and (05i) because one of the 
quadratic forms is singular, independently of the 
Lagrangian. Indeed, for a generic antisymetric 
tensor Hab living in (1 -|- 3) dimensions and an 
arbitrary function /, the following holds 

with 

U = l- - ^{Hab 

1 lo 

* 

and H°‘^= the dual of Note, 

however, that, if Hab is the pullback of a vol¬ 
ume form, we have 

Hab ^t^^'^iLacLbd - LadLbc) = 0, 

meaning that the last term in U always vanishes 
for the integrable models, independently of the 
background solution. Moreover, recalling that 
H^’^Hah = 2 iT 2 , one obtains 

which implies that the characteristic surfaces 
are not given by a product of Lorentzian cones. 
Rather, one of them is given in terms of the zeros 
of a singular quadratic form in Thus, 

we conclude that, for all possible Lagrangians 
depending only on 02 , the corresponding Euler- 
Lagrange equations are not of hyperbolic type. 
Therefore, for these nonlinear field theories, in- 
tegrability is not compatible with hyperbolicity. 
As a consequence, in spite of the fact that the 


models on two-dimensional targets admit infinite 
sets of exact analytical solutions, their linearized 
versions are always problematic. In other words, 
one does not expect that arbitrary initial data 
will launch a nice solution, even locally in time. 

V. CONCLUDING REMARKS 

Let us conclude by mentioning a number of 
possible issues, future directions and open ques¬ 
tions. First of all we stress that, for the inte¬ 
grable models taking values on a bi-dimensional 
target, the associated causal structure is always 
ill-behaved. In other words, the characteris¬ 
tic surfaces do not have the topology of convex 
cones, as expected for a hyperbolic theory. This 
is potentially worrying since small disturbances 
about any (possibly analytic) smooth solution 
will not propagate in a well behaved manner in 
spacetime. What is more, our results do not 
depend on the specihc choice of the Lagrangian 
neither on the topology/geometry of the target 
space. Rather, they indicate that the ve^ notion 
of integrability in the lines of Q] and is not 
compatible with hyperbolicity (well-posedness), 
at least for a 2-dimensional target space. Physi¬ 
cally, this result is somehow expected since these 
theories are effective and often need corrections. 

One may wonder whether slight modifica¬ 
tions of the theory could lead to a hyperbolic 
theory, thus entailing a well-posed Cauchy prob¬ 
lem. It first comes to mind the addition of a 
potential term in the Lagrangian, Unfor- 


9 


tunatelly, a closer inspection of Eq. ([8]) reveals 
that a potential do not contribute to the prin¬ 
cipal part of the former PDEs, implying that it 
is not able to cure the degeneracy of the char¬ 
acteristic surfaces. Another possibility would be 
to consider general Lagrangians of the form 


£ = ClUl -I- 02(72 - U {(f) 


the sextic term in its Lagrangian has connections 
with the pull-back of a volume form to a three- 
dimensional target space, we expect that similar 


results presented here would equally hold, 
is reinforced by a theorem by Wong (see 


'his 


m, 


(26) 


where ci << C 2 , i.e. the Dirichlet term should 
enter as a rather small addition. Eormally, this 
choice is equivalent to the Faddeev-Niemi model 
with a potential and admits nice characteristics 
for sufficiently slowly varying background solu¬ 
tions 27|. However, this model is not integrable 


and it can be considerably difficult to obtain 
exact solutions. Hence, given the fact the in¬ 
tegrable models are physically compelling and 
mathematically simple, it would be extremely 
illuminating to investigate more closely the in¬ 
terconnections between integrability and hyper- 
bolicity. 

An extension worth of future investigation is 
the hyperbolicity of the BPS Skyrme model P- 
11 [. As is well known, the model is based on 
a Skyrme-type low energy effective action which 
does have a Bogomolny bound and exact Bogo- 
molny solutions. It qualitatively reproduces the 
main features of the liquid droplet model of nu¬ 
clei and provides quite accurate binding energies 
of the most abundant higher nuclei (after tak¬ 
ing into account the semiclassical rotational and 
iso-rotational corrections as well as the Coulomb 
interaction and a small isospin breaking). As 


Theorem 9), which implies that the term ui in 
the Lagrangian always introduce a factor that 
is regularly hyperbolic. This in fact has a sta¬ 
bilizing effect on the hyperbolicity of the held 
theory. The inclusion of this term leads us to 
the so called near-BPS Skyrme model, which has 
a chance to be a correct low energy, solitonic 
model of QCD and nuclear matter. We shall 
analyse the evolutionary aspects of this model 
in a forthcomming communication. As a hnal 
remark, we note that the BPS limit realizes an 
important idealization of nuclear matter where 
i) the classical binding energies are zero ii) the 
energ-momentum tensor characterizes a perfect 


huid 


391]. Indeed, the action of the BPS Skyrme 


model is equivalent to the action of a held theo¬ 
retic description of perfect huids in an Eulerian 
formulation 40-^. Perhaps, the degeneracy of 


the characteristic surfaces is related to the ab¬ 
sence of dissipation in the above model or to the 
violation of some of the energy conditions. Up to 
this moment we have no other hints about this 
unexpected behaviour. It would be extremely 
interesting to have any deeper insight into this 
problem. 
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